Strong-coupling theory of superconductivity in a degenerate Hubbard model 
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In order to discuss superconductivity in orbital degenerate systems, a microscopic Hamiltonian is 
introduced. Based on the degenerate model, a strong-coupling theory of superconductivity is devel- 
oped within the fluctuation exchange (FLEX) approximation where spin and orbital fluctuations, 
spectra of electron, and superconducting gap function are self-consistently determined. Applying the 
FLEX approximation to the orbital degenerate model, it is shown that the d^2_j^2-wave supercon- 
ducting phase is induced by increasing the orbital splitting energy which leads to the development 
and suppression of the spin and orbital fluctuations, respectively. It is proposed that the orbital 
splitting energy is a controlling parameter changing from the paramagnetic to the antiferromagnetic 
phase with the ci3,2_j^2-wave superconducting phase in between. 



I. INTRODUCTION 

Recently, new heavy fermion superconductors CeTIns 
(T=Rh, Ir, and Co) have been discovered. These com- 
pounds have been atracted much attention, since vari- 
eties of ordered states are observed by changing tran- 
sition metal ions. Among them, Celrlns and CeCoIus 
are superconducting at ambient pressure with transition 
temperatures Tc=0.4K and 2.3K, respectively.^'^ In par- 
ticular, CeCoIns shows the highest superconducting tran- 
sition temperature among Ce-based heavy fermion sys- 
tems. On the other hand, CeRhlus exhibits an antiferro- 
magnetic transition at a Neel temperature Tn^S.SK and 
becomes superconducting only under hydrostatic pres- 
sure larger than 15 kbar.'^ 

The important experimetal results of these materi- 
als are summarized as follows. Reflecting the fact 
that CeTIns has the HoCoGas-type tetragonal crystal 
structure, quasi two-dimensional Fermi surfaces have 
been observed in de Haas-van Alphen experiments 
of the compounds, consistent with the band-structure 
calculations."'''^ The specific heat of CeCoIus has shown 
considerably large discontinity at the superconducting 
transition temperature by which the superconductivity 
of this compound is considered to be in the strong- 
coupling regime.^ Concerning the superconducting state, 
nuclear relaxation rate of CeTIns exhibits behavior 
below Tc^'^ and thermal conductivity in CeCoIns is found 
to show a component with four-fold symmetry,® which 
strongly suggest the dx^-y^-wave pairing symmetry of 
the superconducting phase of CeTIns. Furthermore, it 
has been shown that in the alloy system CeRhi_2.Ir:rIn5, 
the superconducting phase appears in the neighborhood 
of the antiferromagnetic phase. ^ These experimental re- 
sults indicate that the CeTIus compounds have similarity 
with high-Tc cuprates. Thus, it is natural to expect that 
the mechanism of superconductivity of CcTIns is similar 
to that of high-Tc cuprates. 

From the band-structure calculations of Cclrlns and 
CeCoIns , one can see that characteristic features of the 



Fermi surfaces of these compounds, such as shape, vol- 
ume, and the 4/-electron weight, are almost the same 
with each other. '*'^'^° On the other hand, it has been 
shown in theoretical studies of high-Tc cuprates that hole 
doping leads to the deformation of the Fermi surface and 
the structure of spin fluctuations around the antiferro- 
magnetic vector is changed. Based on the mechanism 
of dj.2_j,2-wave superconductivity induced by the antifer- 
romagnetic spin fluctuation, variation of Tc in the over- 
doped region of high-Tc cuprates has been explained rea- 
sonably by this scenario. "'^'"'^ It means that hole dop- 
ing is controlling the superconducting transition temper- 
ature through the deformation of the Fermi surface in 
high-Tc cuprates. From the point that the superconduct- 
ing transition temperatures of Celrlns and CeCoIns are 
quite different from each other in spite of the similari- 
ties of the Fermi surfaces, it is an improbable scenario 
that the superconductivity of CeTIns compounds is pri- 
marily controlled by the carrier doping as in the high-Tc 
cuprates. Thus, other controlling parameters for the su- 
perconductivity should be searched in CeTIns even if the 
superconducting mechanism of CeTIns is similar to that 
of high-Tc cuprates. In order to find such a controlling 
parameter of superconductivity in the heavy fermion sys- 
tem CeTIn5, theoretical study from the microscopic point 
of view is highly required. 

However, theoretical studies for the heavy fermion su- 
perconductivity have been almost restricted in the phe- 
nomenological level^'^ because of the following problems: 
(1) It is difficult to treat the dual nature of /-electrons, 
coexistence of both localized and itinerant character, 
in contrast with the d-electron systems where the itin- 
erant picture is a good starting point. (2) Compli- 
cated /-electronic states are formed by the combined 
effect of crystal structure and orbital degeneracy of /- 
electrons, which leads to multiple Fermi surfaces. For 
the first difficulty, it may be a reasonable assumption 
that even strongly correlated states with such a duality 
of /-electrons are still adiabatically continued from the 
state in weakly correlated systems which is the key as- 
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sumption of the Fermi-liquid theory. With respect to the 
second difficulty, it may be appropriate to introduce a 
microscopic model considering the essential part of the 
complicated crystal structure. Furthermore, it may be 
necessary to incorporate the orbital degree of freedom 
of /-electrons into the model, since quasi-particle states 
are reflected by kinds of orbitals. Thus, in order to un- 
derstand superconductivity in the heavy fermion system 
from a microscopic point of view, we should develop a 
microscopic theory based on the Fermi-liquid type the- 
ory using an orbital degenerate model. Especially, in 
view of the large discontinuity of the specific heat at Tc 
in CeCoIns, it is necessary to develop a strong-coupling 
theory for superconductivity in CeTIns compounds. 

In this paper, we focus on the effects of orbital de- 
grees of freedom on superconductivity based on a mi- 
croscopic theory applied to a microscopic model with 
the orbital degeneracy. In the next section, we intro- 
duce the orbital degenerate model obtained by includ- 
ing important characters of CeTIns. Then, in order to 
study the superconducting transition in the orbital de- 
generate model, we develop a strong-coupling theory us- 
ing the fluctuation exchange (FLEX) approximation'^'' in 
which spin and orbital fluctuations, the single-particle 
spectrum, and superconducting gap function are deter- 
mined self-consistently. Finally, we discuss experimental 
results for CeTIns in the light of the present theory. 

II. MODEL HAMILTONIAN 

First let us introduce local basis for /-electron systems. 
As is well known, 14-fold degenerate /-electronic states 
split to j—5/2 and 7/2 multiplets due to strong spin-orbit 
coupling where j is total angular momentum. It is quite 
natural to consider only the lower j=5/2 multiplet con- 
tributes to low-energy excitations. This multiplet splits 
into r7 doublet and Fg quartet under cubic crystalline 
electric field (CEF), and the corresponding eigen-states 
are given by 

|rr.> = V^l4)-/f|T^), (1) 

|rfi) = l±i), (3) 

where \jz) are basis of 7=5/2 multiplet, and +(— ) in 
the subscripts denotes "pseudo-spin" up(down) in each 
Kramers doublet. Furthermore, under tetragonal CEF, 
two Fy and one Fg Kramers doublets are formed. Exper- 
imental data of magnetic susceptibility of CeTIns anal- 
ysed by using the CEF theory seem to be consistent with 

the level scheme where the two F7 are lower than the 
r 15,16 

i 6- 



Here we consider that superconductivity in systems 
with the orbital degrees of freedom is affected primar- 
ily by splitting enegry between lowest and excited states, 
while kind of excited Kramers doublet may play sec- 
ondary role to determine details of electronic properties. 
Based on this belief, in the following, we consider only 
iFg) states. In other words, one F7 state is assumed 

to be the highest energy state. Note that |Fg2) and 
(2) 

iFgji) belong to r7 and Fg irreducible representations, 
respectively, in the tetragonal system. Although this as- 
sumption is not exactly the same as the level scheme 
obtained from experimental results mentioned above, us- 
ing Fg and F7 states instead of two F7 states will pro- 
vide even reasonable result to discuss superconductivity 
of CeTIn5.'^ We stress that the Hamiltonian constructed 
from Fg quartet is the simplest possible model including 
essential physics of interplay between pseudo-spin and 
orbital degrees of freedom. 

We include itinerant features of /-electrons by con- 
sidering nearest-neighbor hopping of /-electrons. Since 
realistic effective hopping of /-electrons includes that 
through hybridization with conduction electrons, such ef- 
fects may be renormalized to the nearest-neighbor hop- 
ping of /-electrons after the conduction electron degrees 
of freedom are integrated out. In the present case, the 
matrix elements of the hoppings depend on not only the 
hopping directions but also kinds of orbitals, since the 

forms of the wave functions of |Fg2) and |Fg2) states 
are different from each other. Noting that CeTIns has 
a tetragonal crystal structure and quasi two-dimensional 
Fermi surfaces have been experimentally observed,^'^ it 
is natural to consider the two-dimensional square lat- 
tice composed of Ce^"*" ions. Considering these points, 
we have estimated the hopping matrix elements through 
the cr-bond (//cr) using the tight-binding method. '^'^'-'^ 
The matrix elements of nearest-neighbor hopping of /- 
electrons t^^, between r and r' orbitals along the a- 
direction are given by 

til = -V3i^2 = -^t^i = 3tj2 = t, (4) 
for a=x, and 

<n = v^i?2 = ^tli = 3<^2 = t, (5) 

for a=y.'^ We use t=l as energy unit in the following. 
By further adding the on-site Coulomb interaction 

terms among /-electrons, an effective Hamiltonian of 
CeTIns compounds with orbital degrees of freedom is 
obtained as 

J2 ^rx'/L/i+aTV-A^(nii<.-ni2.)/2 

iarr'fT i 

+ U'^ nir^riiri + U''^ nii„ni2a' , (6) 

ir 1(7(7' 

where fira is the annihilation operator for an /-electron 

(t) 

with pseudo-spin a in the r-orbital state Fg at site i. 
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a is the vector connecting nearest-neighbor sites, and 

niTa=fiTa-fi'r(7- The first term represents the nearest- 
neighbor hopping of /-electrons. The second term ex- 
presses the tetragonal CEF, represented by an energy 
splitting A between the two orbitals. In the third 
and fourth terms, U and U' are the intra- and inter- 
orbital Coulomb interactions, respectively. We ignore the 
Hund's rule coupling, since it may be irrelevant for the 
quarter-filling case with one /-electron per site. To kcx^p 
the rotational invariance in the orbital space for the in- 
teraction part of the Hamiltonian, U' should be equal to 
U when we ignore the Hmid's rule coupling. Thus, in 
this paper, we restrict ourselves to the case of U=U' . 

Considering property of the hopping matrix element, 
the present Hamiltonian in the quarter-filling may be re- 
garded as an effective model in the hole-picture for the 
Cu02 plane of a parent compound of high-Tc cuprate 
La2Cu04, although the practical value of A may be con- 
siderably large for the c?-electron system. This fact means 
that some of the results obtained by using the present 
Hamiltonian can be used for the cuprate with a suit- 
able choice of the parameter set. We also note that in 
the quarter-filling case, the present model is reduced to 
a half-filled single-orbital Hubbard model in the limit of 
A=oo. 

Here we briefly discuss symmetry properties of the 
present Hamiltonian. Since all pseudo-spin operators 
= (1/2) Ei,. E.,.' fLr^a'hr.' {a=x, y, z) with f« 
being the Pauli matrices commute with this Hamiltonian, 
the system is invariant with respect to the rotation in the 
pseudo-spin space. Since the Hamiltonian is consisted of 
the pairs of creation and annihilation operators corre- 
sponding to each Kramers doublet, it is easily confirmed 
that the Hamiltonian is invariant under the time reversal. 
Since we consider the system of two-dimensional square 
lattice composed of Ce^^-ion, the Hamiltonian commutes 
with all elements of the D4/1 point group. Obviously, the 
Hamiltonian is [/(l)-gauge invariant. Thus, the Hamil- 
tonian has the symmetry of D/^h x SU{2) x 6 x U{1) 
after all where SU{2) describes the rotation group in the 
pseudo-spin space and G the time reversal symmetry. In 
the following, we call pseudo-spin as "spin" for simplicity. 

III. FORMULATION 

In our previous work, we have developed a weak- 
coupling theory for superconductivity based on the same 
orbital degenerate model described above, using the 
static spin and orbital fluctuations obtained within the 
random phase approximation (RPA).^" On the other 
hand, considerably large discontinuity of the specific heat 
at the superconducting transition temperature has been 
observed in CeCoIus.^ Since the discontinuity is much 
larger than the specific heat just above the superconduct- 
ing transition temperature, the mass enhancement due to 
the strong interaction between /-electrons may not be re- 



sponsible for the large discontinuity. Rather, this exper- 
imental fact seems to require the strong-coupling theory 
in order to understand the superconductivity in CeTIns 
compounds. Thus, we should develop a strong-coupling 
theory for superconductivity based on the orbital degen- 
erate model. 

In the present paper, we apply the fiuctuation ex- 
change (FLEX) approximation^* to the orbital degen- 
erate model discussed in the preceding section. The 
FLEX approximation provides the Dyson- Gorkov equa- 
tion where the normal and anomalous self-energies are 
obtained on an equal footing, namely a kind of the strong- 
coupling theory. Here, it should be noted that the strong- 
coupling theory includes two major changes for the or- 
bital degenerate system compared with the single orbital 
case. One is multi-component nature of the supercon- 
ducting order parameters in the orbital space, namely or- 
bital symmetric and orbital antisymmetric gap functions. 
The other is the effect of mode-mode coupling among 
spin and orbital fluctuations. We expect that the mode- 
mode coupling significantly affects the temperature and 
frequency dependences of spin and orbital fluctuations. 
Therefore, in order to discuss superconductivity induced 
by these fluctuations the mode-mode coupling effect is 
very important. In the following, first we discuss general 
relations satisfied by the Green's functions required from 
the symmetry of the system, then develop the scheme of 
the FLEX approximation for the degenerate model. 

A. Definition and Properties of Green's Function 

The normal Green's functions G^„(k, t) describing the 

propagating process of electrons from n-orbital to m- 
orbital with moment k and spin a and the anomalous 

ones F^'^ (k, t) and F'!^^ (k, r) describing the supercon- 
ducting condensation are defined as 

G^Jk,T) = -(T.[/w(t)/L,(0)]), (7) 
F^^' (k, r) = -{Tr[hmAr)f-Una' (0)]), (8) 

fZAKt) = -(T.[/lk_(r)/t_,(0)]), (9) 

where hma{T) = e(^-''^)^/k„.e-(^-^^)^ with r being 

the imaginary time, N = J2i,rn,a flnafi^'^ operator 
of the /-electron number, /i the chemical potential, and 
Tt describes the time ordered product. In these equa- 
tions, (• • •) means the thermodynamical average. It is 
convenient to transform the imaginary time Green's func- 
tions to the frequency representation given by 

a(k,ia;;) = / ^ rfTe''^'^a(k,r), (10) 

JO 

where 5(k, r) represents a component of the normal or 
anomalous Green's functions defined above, and u)i={2l+ 
l)7rT is the Matsubara frequency for fermions. 
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Then, under the assumption that the superconducting 
transition does not break the time reversal symmetry for 
the orbital degenerate system, Green's functions satisfy 
the relations as 

GL,(k,ic.O - GZ^iK-iuji)* = GZ^{--k,iuJi), (11) 

= aa'FZ^{-k, -iuji)* = aa'F^n^i-k, iui), (12) 

where the last equalities in these relations arc obtained by 
the time reversal invariancc. Due to the SU (2)-symmetry 
in the spin space, (k, iwj) are decomposed into the 
spin-singlet and spin-triplet irreducible representations, 
defined as 

F^„(k, ioji) = liFiiiK hJi) - FiiO^, icoi)), (13) 

and 

Fl„(k,ia>0 - ^iFliik,iu;i) + Fiiik,iLOi)), (14) 

where F^„(k, iw;) defined above is a representative of the 
three components of the spin-triplet pairs. In the para- 
magnetic system, since we can suppress a superscript de- 
scribing spin state of G^„(k, iw;), we obtain the relations 
for the normal Green's functions as 

Gm„(k,iu;;) = G„m(k, -iwj)* = G„m(k,iwi). (15) 

Furthermore, when the orbitals are defined at the Bravais 
lattice, the inversion operation changes only the wave vec- 
tor k to -k independent on the orbital state as well as spin 
of the quasi-particles forming the Cooper pair. Thus, we 
obtain the relations for F^^{]s., iui) (^=s or t) as 



= F4„(k, -io;;)* = F^mi^i' 



(16) 



Among i4ii(k,iwi), i^„„(k,iwj), and i4„(k,iwi)> we 
also obtain 

F4„(k, iuji) = FT,i(k, iuji) = FiL(k, iuJi). (17) 

Note that these relations are obtained in the time reversal 

invariant system regardless of the spin-singlet and even- 
parity pair or the spin-triplet and odd-parity pair. We 
emphasize that these relations are quite useful to make 
the formulation simple. 

By transforming _F4„(k, iij;) (^=s or t) to the 
"orbital-symmetric" and "orbital-antisymmetric" repre- 
sentations, four real anomalous Green's functions are de- 
fined as 



F«(k,ic^,) = ^ii(k,iu;z), 
F|(k,ia;0 = ^22(k,iwO, 



(18) 
(19) 



F|(k, icji) ^ ^{F^iK i'^i) + Fi.ik, iu;i)), (20) 
F|(k, iui) = -^{F^2{K M) - F|i(k, iwi)), (21) 



where f|(k, iuji) just describes the orbital-antisymmetric 
pairing. For the spin-singlet state, they satisfy the rela- 
tions as 



F^iKiuJi) = F^{-k,icoi) = F^(k,-ia;0, 
F|(k,ia;0 = Fl{--k,iuJi) = -F|(k,-ia;0, 

and for the spin-triplet state, we obtain 

F^(k,ia;0 = -F^(-k,ia;0 = F'^(k,-iuJi), 



FKKiu>i) 



-Fl{--k,iwi) = -Fl{k,-iui), 



(22) 
(23) 



(24) 
(25) 



where rn represents the number of the orbital-symmetric 
component, namely m=l, 2, or 3. One can see that the 
orbital-antisymmetric anomalous Green's function has 
odd-frequency dependence, irrespective of the spin state 
of the Cooper-pair. Thus, due to the odd-frequency prop- 
erty, the orbital-antisymmetric component of the Cooper 
pair may not provide essential contribution to supercon- 
ductivity. In addition to this feature in the frequency 
space, we should also pay attention to the k-dependences 
of F^^{k,iui) because the lattice system and the local 
wave functions should be rotated simultaneously. For ex- 
ample, a symmetry operation C4 of the tetragonal point 
group rotates the wave vector k, and also changes the 
sign of only the wave function of /-electron. The lat- 
ter effect of the symmetry operation leads to the differ- 
ence between k-dependences of the orbital-diagonal and 
the orbital-offdiagonal anomalous Green's functions, in 
order to preserve the symmetry of superconductivity. Af- 
ter all, when the superconducting state belongs to F irre- 
ducible representation, the k-dependences of FfjCkii'^;) 
and F21 (k, ia;; ) have the F x Big-symmetry while the sym- 
metry properties of ^^^(k, iw;) and F22{k.,iuji) in the k- 
space behave as the F irreducible represenation, the same 
as the symmetry of the order parameter. These relations 
obtained in this subsection are useful for practical calcu- 
lations. 

Here we comment on difference between the strong- 
and weak-coupling theories for the degenerate model. 
As will be obtained below, the superconducting gap 
functions are not independent for the orbital-symmetric 
and orbital-antisymmetric parts generally. On the other 
hand, in the weak-coupling theory of the same de- 
generate model, the orbital-symmetric components and 
the orbital-antisymmetric one of Timn (k) define sepa- 
rate superconducting states from each other. Consid- 
ering the result in the present subsection, the complete 
separation within the weak-coupling theory is due to 
ignoring the odd-frequency dependence of the orbital- 
antisymmetric component of the anomalous Green's 
functions F^{]s.,iivi). 



B. FLEX Approximation for the Degenerate Model 

In order to develop a strong-coupling theory of super- 
conductivity in the degenerate system, it is convenient to 



4 



introduce the Dyson-Gorkov equations, which is in the 
matrix form because of the orbital degree of freedom, 
given by 

G{k) = (fc) + G(°) {k) E(^) (fc) G{k) 

+G^°\k)t^^'^\k) Fi{k)\ (26) 
F^{k) = G^°\k) t^^\k) F^{k) 

-G^°\k) t^^'^\k) G{-k)\ (27) 

where Ti^^^k) is a matrix of the normal self-energies, 
S^(^) (fc) describes a matrix of the anomalous self-energies 
for the spin-^ state pairing and the abbreviation k = 
(kjiui) is used. &°\k) is the matrix of the noninter- 
acting Green's function whose matrix elements are given 

by 



GW(k,ia;0 = E 



(28) 



with 



Ekll + ek22 



■P 



Qikll = Q!k22 = 
Q!kl2 = -Q!k21 



1 + 



+ (-1) 

-Ek2 — Ekll 



^, / (£kii -ek22)' 1 ,2 



Ek2 — ek22 
ekl2 



4 

-1/2 



-Ek2 — ek22 



Q!kll, 



and 



Ckii = 2 i (cos kx + cos ky) - — - jj,, 

2 / , , ^ A 

ek22 = 2 i (cos fca; + cos fcj^j + y - ^, 
2 

eki2 = ek2i = --^ t (cosfcx - cosAiy), 



ki2- (29) 
(30) 
(31) 

(32) 

(33) 
(34) 



where £^kp is the energy dispersion of p-th band and akpm 
is the weight of m-th orbital for the p-th band at wave 
vector k. 

In order to obtain a concrete form of the matrix el- 
ements of the self-energy for the degenerate model, we 
adopt the FLEX approximation. First, we define a 
"Luttinger-Ward" functional (A;)] consisting of ladder 
and bubble diagrams for particle-hole processes but ig- 
noring particle-particle processes, and then generate ma- 
trix elements of the self-energy by functional differenti- 
ation of $[5(A;)] with respect to Q{k) where Q{k) rep- 
resents a component of the normal and/or anomalous 
Green's functions. By this procedure, the self-energies 
are obtained within the FLEX approximation by 



S^^LW = ^Y.Y.'^:^,^n{ci)G^.{k-q), (35) 

q nv 



and 



S«l^Hfc) = j^Y.Y.TUnM)FU^-Ql (36) 

where an abbreviation q = (q, i^l;) is used with the bo- 
son Matsubara frequency 0;=2i7rT. The fluctuation ex- 
change interactions used in Stow (A;) are given by 



neflf 



+ U'')f\q){W + U°) + 3U' - UXm,.n. (37) 

For the spin-singlet channel, matrix elements of the ef- 
fective pairing interaction are given by 



+ U°) + W + u°\ 



(38) 



and for the spin-triplet channel, we obtain 
1 



+ \iU' + U°)x^{q){U'' + 



U°)+U' + U°]f,m,n., (39) 



with 



[/ 

[/ 

[/' 

[/' 



U 2U' 
2U' U 









-U' 

-U' 



As has alredy been pointed out in the previous work,^° 
one can see that the contributions to the pairing interac- 
tion of the spin and orbital fluctuations arc, in general, 
destructive for the spin-singlet channel while constructive 
for the spin- triplet channel. X^{q) and X°{q) in the above 
expressions of self-energies correspond to the spin and 
orbital fluctuations, respectively, and within the FLEX 
approximation these are given as 



with 



nQ) = [^-U'xiQ)]-'xiq), 
r{q) = [i + U''x{q)]-'x{q), 

x{(i) = x"'{(i)-x"'{(i), 

^ O - ^ (7(7 , ^ (7(7 f , 

X (?) = X (?) +x (9), 



(40) 
(41) 

(42) 
(43) 



where x (<?) is the matrix of the irreducible susceptibil- 
ity corresponding to one bubble Feynman diagram with 
using the renormalized Green's functions, given by 



X (9) 

XllxM X?f22('Z) X?fl2('?) Kf2l(9) 

^211 (9) X22"22(9) X2212('?) X22"2l(9) 

XlilM X?222(<?) Xm'2(«) X?22l(9) 

Xmi(9) X2f22(?) X2n2(9) X2f2l(9) 



(44) 
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Here we define 



IV. CALCULATED RESULTS 



and 



T 

Wo 



k 



G,,{k + q)Gjt{k), 



T 

iVo 



(45) 



(46) 



Although these expressions are derived for the doubly 
degenerate system, it is straightforward to extend this 
formalism to systems with more orbital degrees of free- 
dom. 

To calculate Tc we linearize the Dyson- Gorkov equa- 
tions with respect to F^{k) or I]^(^)(fc). The transi- 
tion temperature for superconductivity is determined as 
the temperature below which the linearized equation for 
S^(^^(fc) has a nontrivial solution. The linearized gap 
equation is given by 



N, 



q) (47) 



with 



wL,Al)=T.^^^-^,nMG^sik-q)G,tiq-k). (48) 



Finally, we mention actual calculations in the FLEX 
approximation. The FLEX calculation is numerically 
carried out at fixed parameter values of U =U'=4: and one 
/-electron per site on the average. Summations involved 
in the above self-consistent equations are performed us- 
ing the fast Fourier transformation algorithm both for 
the k-space with 32 x 32 meshes in the first Brillouin 
zone and for Matsubara frequency sum. In particular, 
with respect to the Matsubara frequency sum, we have 
adopted an useful method developed by Deisz et al?^ to 
include high frequency contribution efficiently. When the 
relative error of every matrix element of S'^^^(fc) for all k 
and uji becomes smaller than 10~^, we assume that the 
solution is obtained for the self-consistent equations. In 
the present calculations, a second-order magnetic transi- 
tion is defined by 



det[i - Wt (?)] = V. 



(49) 



where -q = 0.002 is used throughout this paper. In- 
troduction of such a small value for rj may be under- 
stood as the effect of weak three-dimensionality which 
is ignored in the present treatment. In general, increas- 
ing three-dimensionality extends the antiferromagnetic 
phase, while it suppresses the superconducting phase. 
A choice of sufficiently small value of rj is considered to 
be consistent with a description of the quasi-two dimen- 
sional system such as CeTIus compounds. 



A. Spin and Orbital Fluctuations 

In this section, we show the results obtained within 
the FLEX approximation described in the previous sec- 
tion. We start from properties of the spin and orbital 
fluctuations of the orbital degenerate model. Here, in or- 
der to characterize the strength of the spin and orbital 
fluctuations, we define and a", respectively, as 



a" = Min det[l-Wx (q,0)], 
a° = Min det[i + C/°f°(q, 0)], 



(50) 
(51) 



where "Min" in these expressions means the minimum 
value in the q-space. These quantities indicate inverses 
of enhancement factors of dominant spin and orbital fiuc- 
tuations in the momentum space. Note that decreases of 
and a° correspond to developments of the spin and 
orbital fiuctuations, respectively, at a wave vector in the 
q-space. In the present approximation, vanishing of an 
eigenvalue of [1 — tj^x (<3i 0)] or [1 + U°x (q, 0)] defines 
an instability of the spin or orbital density fiuctuation. 
Thus, and o° may be used as indicators of the strength 
of spin and orbital fluctuations. The temperature depen- 
dence of and a° are shown in Fig. 1 for various val- 
ues of the orbital splitting energy A. With decreasing 
the temperature, decreases while q° is almost inde- 
pendent of the temperature. By increasing the orbital 
splitting energy A, is considerably suppressed at suf- 
ficiently low temperatures while a° increases. 




0.80 



0.0 0.1 0.2 0.3 0.4 

T 

FIG. 1. Temperature dependences of (a) and (b) a° for 
each value of A. 

In order to investigate the dynamical properties of the 
spin and orbital fiuctuations, we define local spin and 
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orbital susceptibilities as 

^^xtUuH = ]^ E I'^X^.llJq, ^ + (52) 
° q 

where (5 is a positive infinitesimal quantity. The ana- 
lytical continuation of x^°^mn{q) to the real axis is car- 
ried out by the Pade approximants. In order to clarify 
contributions to the longitudinal and transverse compo- 
nents of the orbital fluctuation from Im^J^j^ we 
introduce operators of the charge density, longitudinal 
orbital density, and transverse orbital density as pi =^ 

(1/2) X]m,cr /imcr/i""^' '''i ~ / ^m,n,cr fima'''nm finer , 

andr?" = J2m,,i,a fLai^''+^T'^)mnfina, respectively. Dy- 
namical susceptibilities for these operators are defined as 

/•OO 

X°cH = i/ dte"([pi(<),pi(0)]), (53) 
Jo 

xlH = i dte'-'-'H[Tt{t),rtm), (54) 

^0 

X^(c.) = i/ dte'-'-'H[rht),rfm), (55) 
Jo 

where these quantities correspond to the local compo- 
nents of the net charge fluctuation, the longitudinal or- 
bital fluctuation, and the transverse orbital fluctuation, 
respectively. By using X^iy.mn(^); these dynamical sus- 
ceptibilities are described as 

Xc(^) = liXiiM^) + X22,22(^) + 2Xn,22(^)], (56) 

Xli^) - \[Xii,n{^) + X22,22(^) - 2Xn,22(^)], (57) 
Xt(^) = X21,2i(^), (58) 

where we make use of relation Xii22('^) — X22ii(^)- 
The frequency dependence of Imx^^m„(w) (Imxx('^)) 
for A/t = and A/t = 2 are shown in Fig. 2(a) and 
Fig. 2(b) (in Fig. 2(c) and Fig. 2(d)), respectively. From 
these figures, one can see that the intensity of Imxii^ii(ij-') 
in low-energy region develops with increasing A while in- 
tensities of two components of the orbital fluctuations are 
suppressed and net charge fluctuation almost unchanged. 
In the present case, it is understood that the spin fluc- 
tuation of the /-electron belonging to r=l orbital (cor- 
responding to Fg^'') provides predominant contribution 
compared with other fluctuations. The longitudinal or- 
bital fluctuation seems to play the secondary role. 

Here we comment on the reason why intensities of the 
orbital fluctuations in the low-energy region do not de- 
velop significantly even for the case of A=0. We note 
that any spectral functions satisfy the sum rule in the fre- 
quency space. It means that suppression of the spectrum 
in the low-energy region leads to enhancement of that 
in the high-energy region or vice versa. In the present 
Hamiltonian, the /-electron number in the t=1 orbital 
ni is larger than n2 due to the difference of the strength 



of hopping integrals included in the kinetic term ekn and 
£k22- Since ni — 712 ^ in this case, the orbital polar- 
ization represented by is non-zero, the longitudinal 
orbital fiuctuation should be suppressed by the factor of 
(t/')^. This is the reason why the spectrum of the lon- 
gitudinal orbital fiuctuation in the low-energy region is 
suppressed even for the case of A=0. 

In order to discuss the nature of transverse orbital fiuc- 
tuation, we derive the asymptotic forms of Xm °,mn (q, z) 
for z 3> 1 as 

x:.l°L„(q, z) « -^x^l°L„(q)^ + o(^), (59) 

with coefficients 

^Xi2°l2(q) = -'^X2i°2i(q) = ^("1 - '^2), (60) 

for (^i/,r7in) = (12,12) or (21,21) and 

Sxtlni^i) = (61) 

for other cases. Note that X2i2i(qi^) J^^^ equal to 
X^{ci,z). From these expressions, we can conclude that 
in the present case of ni — n2 ^ 0, the intensity of spec- 
trum for the transverse orbital fiuctuation shifts to higher 
energy region to satisfy the sum rule. Thus, the signif- 
icant difference of the /-electron numbers for the two 
orbitals prevents essential development of the low-energy 
excitation for the orbital fluctuations. By noting that 
the value of ni — 712 is controlled by the orbital splitting 
energy A included in ekii and ek22; it is considered that 
not only the longitudinal orbital fluctuation but also the 
transverse orbital fluctuation are suppressed by increas- 
ing the orbital splitting energy. 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



(D (0 

FIG. 2. Spectra of the local components of (a) spin and 
(c) orbital susceptibilities for A=0. (b) and (d) are for A=2. 
For (a) and (b), solid, dotted, dashed, and dash-dotted lines 
describe components of ni/mn— 1111, 2222, 2121, and 2112, 
respectively. For (c) and (d), solid, dashed, and dotted lines 
correspond to X=L, T, and C, respectively. 
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The momentum, q, dependences of the principal com- 
ponents of x''(q, 0) and x°{c[,0) are shown in Fig. 3, at 
a fixed temperature r=0.02 for different orbital split- 
ting energy, A. The upper and lower panels indicate 
the results for A=0 and A— 2, respectively. For A=0, 
the antifcrromagnetic spin fluctuation in the t=1 orbital 
corresponding to Fg^^ is enhanced, but not sufficiently de- 
veloped to induce d^2_y2--weLve superconductivity. With 
increasing the orbital splitting energy to A— 2, the anti- 
ferromagnetic spin fluctuation for the t=1 orbital fur- 
ther develops, and orbital fluctuations arc completely 
suppressed compared with the well developed antifcrro- 
magnetic spin fluctuation. 
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FIG. 3. q-dependences of (a) spin and (b) orbital suscep- 
tibilities for A=0. (c) and (d) are for A=2. For all figures, 
solid, dotted, dashed, and dash-dotted lines describe compo- 
nents of nfmn^llll, 2222, 2121, and 2112, respectively. 

Let us compare the present results shown in Fig. 3 with 
those obtained within the RPA in the previous work.'^'^ 
For A=0, the antifcrromagnetic spin fluctuation in the 
T=l orbital already develops without increasing A, while 
similar momentum dependences for many components of 
the spin and orbital susceptibilities has been seen within 
the RPA. When the orbital splitting energy is increased, 
the antifcrromagnetic spin fluctuation in the t=1 orbital 
obtained within the FLEX approximation is considerably 
enhanced in comparison with that within the RPA. The 
considerable development of the antifcrromagnetic spin 
fluctuation may rather prevent dx^-y^-^&ve supercon- 
ductivity because of decoupling effect through the fluc- 
tuation exchange self-energy. 



B. Phase Diagram 

In Fig. 4, the phase diagram obtained within the 
FLEX approximation for the orbital degenerate system 
is shown, where the solid and open circles describe the 



superconducting and antifcrromagnetic transition points, 
respectively. Although we do not carry out the FLEX cal- 
culation in the broken symmetry states, we may expect 
that the first-order phase transition between the d^^^y^- 
wave superconducting phase and the antifcrromagnetic 
phase takes place in the quasi-two dimensional system 
as discussed for the three-dimensional Hubbard model. 
For the latter case, it has been shown that the increase 
of the three-dimensionality expands (shrinks) the mag- 
netic (ci2,2_j,2-wave superconducting) phase and the co- 
existent phase between the d^2_y-2-^ave superconducting 
and magnetic phases appears only for the system having 
the moderate three-dimensionality. From this considera- 
tion, the dotted curve is drawn as the first-order phase 
boundary expected between the two phases. 
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FIG. 4. Phase diagram in the T-A plane for !7=!7'=4.0 
obtained by the FLEX approximation. The dotted curve is 
drawn by hand as expected phase boundary of first-order 
phase transition between the d^2_j,2-wave superconducting 
phase and the antiferromagnetic phase. In the figure, PM, 
AF, and SC(Big) express paramagnetic, antiferromagnetic, 
and d^2_y2-wave superconducting phases, respectively. 

From Fig. 4, we see that (1) the spin-singlet supercon- 
ducting phase with Big-symmetry appears next to the 
antiferromagnetic phase and (2) Tc is enhanced with in- 
creasing the orbital splitting energy A. Recalling the 
expression of the effective pairing interaction for spin- 
singlet state T'^mm^(9)j the increase of Tc with increas- 
ing A is ascribed to the development of the antiferro- 
magnetic spin fluctuation and the suppression of the or- 
bital fiuctuation, where the spin and orbital fiuctuations 
are destructive each other for the spin-singlet pair. From 
these observations, we conclude that the superconducting 
phase is induced by the development of the antiferromag- 
netic spin fiuctuation for the r=l orbital accompanied by 
suppression of the orbital fluctuation with increasing the 
orbital splitting energy A. 
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V. DISCUSION AND SUMMARY 

Based on the simple model with orbital degree of free- 
dom, we have shown that in the system of the density of 
one electron per site, the (ia;2_j^2-wave superconducting 
transition temperature becomes higher for the larger or- 
bital splitting energy. For sufficiently large orbital split- 
ting energy, the antiferromagnetic transition takes place 
at a Necl temperature. Here we discuss the experimen- 
tal results for CeTIns from the present theoretical re- 
sults. Analyses of experimental data of magnetic suscep- 
tibilities by using the CEF theory have determined the 
level schemes of CeTIns compounds where the two Fy 
are lower than the Tq}^'^^ The energy splitting between 
the two F7 is estimated as 68K for CcRhlns, 61K for 
Celrlns, and 151K for CeCoIns. As we have mentioned in 
introducing the present model, we focus on these orbital 
splitting energies of the compounds regardless of kind 
of excited Kramers doublet. Considering that the effect 
of the orbital degree of freedom is generally quenched 
by increasing the orbital splitting energy, the present re- 
sult may be independent of details of level scheme in 
orbital degenerate models. Namely, the present result 
seems to be consistent with the experimental one in the 
sense that Tc in CeCoIns with larger orbital splitting en- 
ergy is higher than that in Celrlns. Thus, one can expect 
that the orbital splitting energy actually plays a role of 
controlling parameter for the superconductivity around 
the antiferromagnetic phase in the heavy fermion system 
of CeTIns compounds. 

On the other hand, the present theory seems to be in- 
adequate to explain the antiferromagnetism of CeRhlns 
with the level scheme similar to Celrlns. Recent de 
Haas-van Alphen experiment of CeRhlns has reported 
that the Fermi surfaces are almost unchanged up to 2.1 
GPa slightly higher than the critical pressure at which 
the curves of the siiperconducting and antiferromagnetic 
boundaries cross each other. ^'^ This may mean that the 
/-electronic states are much lower than the Fermi level, 
and local character of /-electron is dominant in this com- 
pound. Thus, it will be a future problem to clarify the 
antiferromagnetism of CeRhlns by developomg a micro- 
scopic theory based on the local property of /-electrons. 

It is instructive to recall that the present Hamilto- 
nian with considerably large A in the quarter-filling is 
regarded as an effective Hamiltonian in the hole-picture 
for the Cu02 plane of the cuprate La2Cu04 which is 
the parent compound of high-Tc cuprate where d^^-yi- 

and d3j,2_^2-orbitals just correspond to Fg^^ and Fg^^ in 
the present model, respectively. Considering this similar- 
ity, it may be a reasonable speculation that a supercon- 
ducting phase may appear next to the antiferromagnetic 
phase of La2Cu04 in some situation where distance be- 
tween the apical O- and the planer Cu-sites is shrunk, 
the orbital splitting energy between dx2-y2- and d^z^-r^- 
orbitals is decreased, after all. 

In summary, based on the effective microscopic model 



with orbital degeneracy for /-electron systems, we have 
developed the strong-coupling theory for superconductiv- 
ity. Considering the time reversal symmetry involved in 
the model, one can show that the orbital antisymmetric 
component of gap function hardly contributes to super- 
conductivity because of the odd-frequency dependence. 
Using the FLEX approximation in which the spectra of 
/-electrons, the spin and orbital fluctuations, and the su- 
perconducting gap functions are determined consistently, 
it has been shown that the dx2_y2-wa,ve superconducting 
phase is induced by increasing the orbital splitting en- 
ergy which leads to the development and suppression of 
the spin and orbital fluctuations, respectively. Based on 
these results, we have proposed that the orbital splitting 
energy is the key parameter controlling the changes from 
the paramagnetic to the antiferromagnetic phases with 
the rfj,2_y2-wave superconducting phase in between. 
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